Stabilization of Majorana modes in vortices in the superconducting phase of 
topological insulators using topologically trivial bands. 
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It has been shown that doped topological insulators, up to certain level of doping, still preserve 
some topological signatures of the insulating phase such as axionic electromagnetic response and the 
the presence of a Majorana mode in the vortices of a superconducting phase. Multiple topological 
insulators such as HgTe, ScPtBi, and other ternary Heusler compounds have been identified and 
generically feature the presence of a topologically trivial band between the two topological bands. 
In this letter we show that the presence of such a trivial band can stabilize the topological signature 
over a much wider range of doping. Specifically, we calculate the structure of vortex modes in the 
superconducting phase of doped topological insulators a model that capture the features of HgTe 
and the ternary Heusler compounds. We show that due to the hybridization with the trivial band, 
Majorana modes are preserved over a large, extended doping range for p-doping. In addition to 
presenting a viable system where much less fine-tuning is required to observe the Majorana modes, 
our analysis opens a route to study other topological features of doped compounds that cannot be 
modeled using the simple Bi2Se3 Dirac model. 
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Solid-state realizations of Majorana fermions, parti- 
cles which are their own-anti-particles, are much sought 
after for their promise of new fundamental phenomena 
and associated quantum computing applications [1]. By 
now the list of candidate systems to realize these parti- 
cles has grown quite long including: non-Abelian frac- 
tional quantum Hall states 2J, chiral superconductors ^Sj, 
topological insulators in proximity to superconductors [1], 
axion strings [5], and low-dimensional spin-orbit cou- 
pled semiconductors coupled to superconductors [6l [7]. 
One of the most promising directions relies on the in- 
terplay between topological insulators (TIs)[S] and s- 
wave superconductivity [31 [S]- Two such routes are: 
(i)the trapping of Majorana bound states (MBS) in 
vortices found in Tl/s-wave superconductor proxim- 
ity heterostructures[4; (ii)vortex bound states of doped 
TIs which exhibit an intrinsic s-wave superconductor 
phase [S]. The second mechanism has perhaps greater 
potential, as things currently stand, since it does not re- 
quire the fine-tuning of TI material properties to make 
the system a bulk insulator ^. Indeed, there is already 
a candidate material: Cu doped Bi2Se3 which is a super- 
conductor below 3.8 K[Tni[IT] although the precise nature 
of the superconducting state is still undetermined. 

The initial proximity effect proposal of Fu and Kane 
illustrates that Majorana zero-modes will be localized 
where superconducting vortex lines in Tl-superconductor 
heterostructures intersect the topological insulator sur- 
face states[4]. The bound states are confined to the vor- 
tex cores but they penetrate into the nominally gapped 
bulk-region below the surface. The localization length of 
the bound states in the bulk naturally depends on the 
size of the bulk gap and if the bulk is doped then the 
bound states can tunnel away from the surface and be 
destroyed^. While this system, in most respects, is a 
trivial s-wave superconductor, it still remembers it arose 
from a TI parent state as long as the doping is not too 



high. The consequence of this topological signature is 
that vortices will still bind MBSs at the intersection be- 
tween the vortex line and the surface. As the TI becomes 
more doped a vortex quantum phase transition (VQPT) 
occurs at a critical chemical potential and beyond which 
no Majorana modes are bound to the vortices. This 
mechanism comes with its own challenges as one desires 
the formation of an s-wave gap (which may or may not 
be favored), and there is only a finite doping range over 
which the system will exhibit this phenomenon. The lat- 
ter restriction leads to the competition between a desire 
to stabilize a superconducting phase (higher doping) and 
preserving the MBSs (lower-doping) and thus might still 
require fine-tuning of the chemical potential. 

Here we find a more attractive material context for 
this mechanism which removes the fine-tuning constraint 
by considering materials with electronic structure simi- 
lar to bulk HgTe. This requires the consideration of 6- 
bands instead of the 4-band model of the Bi2Se3/Bi2Te3 
family [12]. Our essential insight is that the coupling of 
the Dirac structure to an additional "trivial" band(s) can 
strongly affect the location of the critical chemical poten- 
tial and stabilize the MBSs over a much larger range of 
doping by extending the energy range in which the topo- 
logical signature persists. The set of materials where our 
analysis can be applied is quite large and includes bulk 
HgTe, which was recently confirmed to be a TI[T3], and 
the ternary Heusler materials {e.g. ScPtBi) which were 
recently predicted to be TIs[Tlj. Beyond this application 
our results indicate that the presence of a topologically 
trivial band can aid the persistence of topological features 
{e.g. magneto-electric response) beyond the un-doped 
regime. For these materials this implies the existence of 
a unique metallic phase over a broad doping range that 
has not been fully characterized. 

In this letter we will primarily consider bulk HgTe, 
a zero-gap semiconductor that exhibits a TI phase un- 



der compressive strain although our results directly ap- 
ply to the ternary Heusler compounds with just a change 
of parameters in the model [131 IHl US]- The low-energy 
band structure is time-reversal invariant and consists of 
three doubly-degenerate bands near the F-point. The 
bands closest to the Fermi-level nominally consist of the 
J = 3/2 multiplet of a spin-orbit split p-orbital (denoted 
by Fg) which, when (un)strained is (4)2-fold degener- 
ate at the F-point. The other relevant band nominally 
consists of an s-orbital (Fg) and forms a Dirac Hamil- 
tonian with the Jz = ±1/2 states (the light hole(LH) 
band) of the Fg multiplet [IS]. The relative orientation 
of the Fg and Fg bands are inverted and thus induce 
a topological band structure [16]. The other doubly 
degenerate band, which is made from the J^ = ±3/2 
states (the heavy hole (HH) band) hybridizes with both 
of the Dirac bands, but does not qualitatively affect the 
topological properties [15] ■ The 6-band model Hamilto- 
nian expanded near the F-point in a basis convenient for 
our analysis [|Fg, 3/2), |F6,l/2), |F8,-l/2), |Fg,-3/2), 
|F6,~l/2), |F8,l/2)]is: 
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where we have introduced a model parameter A which 
we will use to tune the coupling between the HH bands 
and the other bands to pass from the Dirac limit (A = 0) 
to the physical limit (A = 1). Additionally, S is the 
gap generated by the compressive strain, and the re- 
maining parameters are determined by the material band 
structure [TTj. For all of the materials in which we 
are interested E^ > and E > 0. The Hamiltonian 
iJugTe preserves time reversal symmetry with the opera- 
tor T — —i(Jy ® 13x3© where Q is complex conjugation. 
When kz = 0, ffngTe is block diagonal with two decou- 
pled 3x3 blocks which are time-reversed partners. The 
energy spectra and orbital composition when A = 0(1) 
are presented in Fig. la(d). 

To determine the existence of MBSs on superconduc- 
tor vortices we will add s-wave superconductivity to 
the 6-band model at the mean-field level and study the 
Bougoliobov-De Gennes (BdG) spectrum when a vor- 
tex line is oriented along the z-direction. Note that we 
are not considering proximity induced superconductivity 



but s-wave superconductivity generated due to a Fermi- 
surface instability in doped topological insulators [9] . Be- 
sides the BdG doubling of the Hamiltonian, one needs 
to include the vortex winding in the phase of the super- 
conducting order parameter. Considering the extreme 
type-II superconductor limit we ignore the magnetic field 
which does not qualitatively change the results [H]. In 
an ordinary s-wave superconductor, the vortex core states 
are gapped [TH] with a "mini" -gap 5 oc ^ where A is the 
superconducting gap and Ep is the Fermi energy of the 
normal metal state. Unexpectedly, it was found that in 
doped topological insulators, these modes can become 
gapless at a critical doping level /Zc leading to the afore- 
mentioned VQPT [S] . Essentially, when tuned away from 
such a critical point, one can think of the vortex line it- 
self as a gapped ID superconductor which can either ex- 
ist in a topological (/ij < /i < /j-lt) or trivial (/x > /j,+ or 
/i < /i^ phase where /ij is the critical chemical poten- 
tial in the conduction and valance bands respectively) 
phase analogous to the Kitaev p-wave superconducting 
wire [3 HO]- When the vortex is in a topological phase it 
will have MBSs localized on its ends i.e. at places where 
the vortex line terminates on a surface. This is exactly 
the phenomenon we seek and thus our goal is to find the 
region of chemical potential where the ID vortex line is 
in a gapped topological phase. Here, instead of studying 
the appearance of Majorana modes on the boundary, we 
study the VQPT in a torus geometry with a vortex line 
along the periodic z-direction and use Lanczos exact di- 
agonalization calculations find the lowest energy vortex 
mode [Idl]. 

For comparison we show two different results, each for 
a different couphng A = 0, 1. In the case A = (Fig.lc, 
dashed line) we essentially reproduce the results of Ref. 
m as this is the Dirac limit. We see two critical points 
one in the conduction band and the other in the valence 
band(s); both critical points occur past the onset of bulk 
doping. For A = 1 (Fig. Ic solid line) the HH band is fully 
coupled and we can see a critical point in the conduction 
band that has been pushed toward the conduction band 
edge while the critical point in the valence band has been 
driven to a much lower energy i.e. vastly increasing the 
range of valence band doping which would contain stable 
MBSs. 

This large effect is rather surprising, but we will now 
argue that it is the natural outcome of coupling to a 
trivial band. From the arguments of Ref. [HI we know 
that in the weak-pairing limit (A << Ep) we can de- 
termine the VQPT from the normal-state Hamiltonian 
alone, without using the BdG formalism. In this case, for 
a vortex-line oriented in the z-direction, the critical chem- 
ical potential is determined solely by the Berry phases of 
the normal-state Fermi-surface (FS) in the fc^ = and 
kz — t: planes. For models of the type we are considering 
there are no low-energy states near /c^ = tt so we will 
only consider the k^ = plane. It was shown that for a 
single (spin degenerate) FS that the lowest vortex core- 
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FIG. 1. Band structure and orbital content of the bands of HgTe (indicated by colors: blue jrg,±l/2), red ]r6,±l/2), green 
|r8,±3/2)) for A = (a) and A = 1 (d). Insets show zoomed out energy spectrum. Corresponding Berry-phase on the Fermi 
surfaces at fcz = as a function of chemical potential for A = (b) and A = 1 (e). The Berry phases are calculated by a 
numerical line-integral around each independent Fermi-surface. The dashed green/solid green lines in (e) represent the two 
different HH Fermi-surfaces. Inset shows zoomed out Berry phase so that the tt Berry phase of the Fe band can be seen, (c) 
Exact diagaonlization results for the lowest energy vortex core states for A = 0, 1 dashed/solid respectively where dots indicate 
critical points. This was calculated for fc^ = for a lattice size of 180 x 180 with HgTe parameters adapted from Ref. 1171 and 
S = 0.1, A = 0.3. (f) Valence-band \i^ as a function of A. 



state energies are proportional to En oc [$b -|- (2n — l)7r], 
where n is an integer, so the vortex-line is critical when 
the Berry-phase $b around the FS is tt. Additionally, 
for a four-band Dirac Haniiltonian, there is an easy way 
to determine the energies for FS's with tt Berry phase 
as the fc-dependent mass term vanishes exactly at that 
chemical potential. That is, at this chemical potential 
there is a VQPT point where a zero energy vortex mode 
emerges. For our 6-band model there are some chemical 
potentials where there are more than one independent FS 
when doped into the valence band(s), but let us first con- 
sider cases where there is only a single (spin degenerate) 
FS. While the the presence of zero energy mode directly 
depends on the the value of $s [9^ even the single-FS 
case does not allow us to determine $b = tt in a simple 
way because the vanishing of the mass term of the Dirac 
block does not guarantee $b — tt when the HH band is 
added (A ^ 0). To determine $5 we need to consider 
kz = where Eq. [I] breaks up into two decoupled 3x3 
blocks. Let us consider one block of the Hamiltonian in 
the kxky-plane in polar coordinates (fcy, 9) in the basis of 
[|mj = 3/2), |mj = l/2), |mj == -1/2)] , 
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Due to time reversal symmetry, the Berry phases of H^ 
and H^ only differ by signs so we only consider H^ 
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and calculate the Berry phases at the FS's. To sim- 
plify the calculation, we have used the spherical approx- 
imation which implies that Ai = A2, Bi ^ B2, and 
F = in Eq. [T] and that in-plane FS's will be circu- 
lar. The generic solution can be obtained by using an 
ansatz which captures the 0-dependence and is of the 
form: ^I^ = (a(A:||)e~*^, 6(fc||), c(fc||)e*^) where the vec- 
tor (a, 6, c) is an eigenvector of il-|_(fc||, 0). Therefore, for 
the case of a single in-plane FS, the Berry phase of the 
wavefunction ^ is 
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This result applies in the entire regime of n-doping, and 
in several regions of p-doping where only a single HH FS 
persists. The expression for ^b is generic for our model 
and the weights a,b,c depend on the particular FS, the 
value of yU, and the value of A. 

Let us now illustrate how turning-on the coupling to 
the the trivial HH band affects the locations of critical 
FS's i.e. places where $b = tt. First, we consider A = 
so that the system is in the Dirac limit (see Fig. lb). 
If the Fermi-level is in the HH band the Berry phase 
is an integer multiple of 2tt (since |ap = 1, |cp = 0) 
and thus does not generate a VQPT. However, if the 
Fermi-level is in one of the two Dirac bands then the 
Berry phase is ^b = 27r|cp which will reach a value of 



TT when |c|2 = |6|2 = 1/2 as found in Ref. 0. This 
will occur at one of the chemical potentials n^ (A = 0) = 



Ao+Bx 



' ± J gf'^^Ig^)' for the upper and lower bands re- 
spectively as indicated by the solid dots in Fig. lb. With 
a small A 7^ there will be a correction to the LH and Fg 
band wavefunctions of the form 5^± — ((5±, 0, 0)^ where 
the ± indicate upper(LH)/lower (Fg) band respectively. 
The resulting wavefunctions and Berry phases at the orig- 
inal critical energies are ^± = ((5±, ±1, l)"^/\/2 -I- |(5±p 
and $B± = 7r(l - 3|5±p/(2 -I- |i5±P)) which have clearly 
been decreased independent of which band. The effect 
on the respective critical potentials in each band is, how- 
ever, asymmetric. We can see this by noting for A = 0, 
$B shifts from 2tt to (0 to 27r) as one raises (lowers) 
the chemical potential through the LH (Fg) band (see 
Fig. lb). Our restriction to a single-FS implies that the 
critical points are determined by $s — -k. Thus, with 
respect to the original critical points we see that the per- 
turbed critical points are at ^^* < ^J. This means that 
the viable doping range has decreased for the conduction 
band, but increased for the valence band. Indeed, we 
see it is exactly the coupling to the trivial band which 
causes the shift in the Berry phase. If we assume axial 
rotation symmetry one can easily extend this result to 
trivial bands with higher angular momentum by replac- 
ing \mj = 3/2) with \mj = (n + 1/2)) which will lead to 
$B = 27r(|cp— n|ap). For n > 0, this will share the same 
qualitative features as the mj ~ 3/2 case; but for n < 0, 
the range of the MBSs will decrease for the valence band 
and increase for the conduction band. For most ranges 
of p-doping, when A 7^ 0, this calculation is merely sug- 
gestive since there will be more than one FS, but we will 
see in the next paragraph that these arguments are still 
approximately valid in most circumstances. We also note 
that in Fig. lc,f we have already seen that the critical 
point can be shifted dramatically in the valence band due 
to the nature of the coupling to the trivial HH band, but 
one must eventually take other bands into account, be- 
cause as ^~ is pushed lower and lower we will eventually 
have to consider the spin-orbit split off band which we 
have ignored. The qualitative trend, however, is quite 
striking and is favorable for experimental realization. 

Finally, to complete our analytic interpretation of the 
numerics we consider the important case when two FS's 
are present simultaneously for some range of chemical 
potentials. For our model there are two FS's present in 
two separate energy ranges (when A 7^ 0): (i) a chemical 
potential near the top of the HH band leading to a nar- 
row energy range with two HH FS's (ii) a wide range of 
chemical potentials with one FS from the HH band and 
one from the Fg band. If momentum is approximately 
conserved in the presence of the vortex i.e. if the vor- 
tex profile is smooth and broad, we can treat the two 
FS approximately independently. For case (i) we can see 
from Fig. le that there are a total of two places where 



$B = ±7r one for each FS. These critical points are quite 
unstable and usually couple together to annihilate. On 
the other hand, for case (ii) we see from Fig. le that 
in this energy range the HH FS will never have a point 
where <^b = tt so the core states generated by this FS will 
always remain gapped. So in this regime, in the decou- 
pled FS limit, the VQPT will continue to be controlled by 
the FS of the Fg band and we can continue to apply our 
previous analytic arguments above. To more accurately 
match our numerics we must consider the hybridization 
between the vortex core modes on the two FS's through 
perturbation theory. If we let Ei {E2) be the lowest en- 
ergy core state from the HH (Fg) FS we can write down 
the effective Hamiltonian 
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where a and /3 are the only couplings allowed by the 
BdG particle-hole symmetry. For the decoupled limit 
a = /? = 0, the lowest core-state energy from each sin- 
gle FS can be determined by the the Berry phase of 
the normal-state FS in the fc^ = plane; for case(ii) 
El ex ($1^ + tt) > and E2 oc (^^'^ + tt), where $|" 
and ^g are the Berry phases of the HH and Fg FS's 
respectively. Therefore, E2 = exactly when /i = fi^ 
so that ^g = — TT. When the inter-FS couplings are on, 
HcS has a pair of zero eigenvalues, i.e. a critical point, 
when E1E2 — \a\'^ — \(3\'^ ■ The values of a,/3 depend on 
the details of the vortex core states and are sensitive to 
finite-size effects. Since E2 decreases from positive to 
negative as /z gets lower then if \a\ > |/3| {\a\ < |/3|) the 
critical chemical potential is modified and driven toward 
(away from) the band edge. In our numerical simulations 
we find the former case where the critical chemical poten- 
tial is reached before a Berry phase ^g — — tt is reached. 
While this effect moves the critical point back toward the 
band edge it still does not counter-act the much larger 
shift due to the hybridization between the Fg and HH 
bands as we have seen in Fig. lc,f. We note in pass- 
ing that we have also considered the effects due to bulk 
inversion asymmetry and find they do not qualitatively 
alter our results. 

Our calculations indicate that materials with Dirac 
bandstructures that hybridize with trivial bands can sup- 
port vortex MBSs over much larger ranges of doping. The 
hybridization with the trivial band delays the VQPT in 
the valence band while accelerating it in the conduction 
band. We are optimistic that these effects can be used 
to find and design an ideal material to support Majorana 
vortex states. More generally we showed that the trivial 
band coupling extends the range of the topological signa- 
ture. This merits further investigations on other features 
of doped TIs which are sensitive to the band inversion of 
the parent insulator state. 
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